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Abstract 

We study the prompt production of the Xc(l + ) and Xb( 1 + ) mesons at high energies. 
Unlike x(0 + ,2 + ) production, x(l + ) mesons cannot be created at LO via the fusion of 
two on-mass-shell gluons, that is gg —> Xc,b(l + ) are not allowed. However, the available 
experimental data show that the cross sections for Xc(l + ) and Xc(2 + ) are comparable. We 
therefore investigate four other x(l + ) production mechanisms: namely, (i) the standard 
NLO process gg —> Xc,b( 1 + ) + <7, (h) via gluon virtuality, (iii) via gluon reggeization and, 
finally, (iv) the possibility to form Xc,b( 1 + ) by the fusion of three gluons, where one extra 
gluon comes from another parton cascade, as in the Double Parton Scattering processes. 


1 Introduction 

It is well known that according to the generalised Landau-Yang selection rule a spin-1 meson 
cannot be produced in the fusion of two identical massless spin-1 particles [I]. Therefore the 
inclusive Xc(l + ) cross section offers the possibility to probe some non-trivial dynamics of the 
NLO interaction. 

There are two ways to overcome the Landau-Yang selection rule: either to account for 
different virtualities of incoming gluons (that is to violate their identity) or to consider the 
formation of x(l + ) by three gluons. The third gluon may be emitted as the new secondary 
particle gg —» X(l + ) + 9 process, or it may occur in the initial state in the fusion process 
9 + ( 99 ) —> X fusion. Recall that, as shown in [2j, the role of such three gluon mechanisms 
of J/i/j -meson production increases with energy])] and may even be dominant at the very high 

1 Unfortunately there was a confusion in [2j - for the 3-gluon mechanism the cross section, and not the 
amplitude, grows as Ins. Thus the cross section at LHC energies will be a few times smaller. 
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energies. Recently the three gluon J/if) production was also studied in [3]. Note that, in 
comparison with the J/-0, the dynamics of y c (l + ) formation is richer since, due to its negative 
C-parity, the J/ijj meson cannot be produced in the fusion of two gluons, even if we account 
for the different virtualities of these gluons. 

The available experimental data on prompt production of y c mesons in high energy hadronic 
collisions (see |4j for a recent review) concern, as a rule, the transverse momentum p t distribution 
of the y-meson production rates or the ratios of cr(y Ci 6(2 + ))/cr(y Ci b(l + )) as a function of p t [^] 
The ratio i? 2 i = y(2 + )/y(l + ) is of the order of one and practically does not depend on p t . 
It can be described by a constant, f? 2 1 , independent of p tl both for y c and Xb processes: in 
particular for \b the value R 21 = 0.85 ± 0.07 [5jj^] This was not expected. The lowest-order 
theory predicts the growth of this ratio for decreasing p t , since, at low p t , y(2) meson can be 
produced at LO while y (1) only occurs at NLO (see, for example, [7]). Moreover, in fixed-target 
7 1 — Be interactions at 515 GeV/c the p t - integrated yields of y c (l + ) (a = 232 ± 37 ± 37 nb) and 
y c (2 + ) (<7 = 407 ± 71 ± 69) nb [8j were measured. This does not indicate the expected strong 
suppression of y c (l) which cannot be produced at LO via on-mass-shell gluon-gluon fusion. It 
is worth mentioning that new important information could come from fixed-target experiments 
using an LHC beam, as advocated in [9], where lower p t values could be reached. Thus it is 
topical to study theoretically the different possible mechanisms of y(l) production. 

The inclusive cross sections were calculated at lowest order in ra. Numerical evaluations 
at the LHC energy are given, for example, in [7j. A specific higher-order process - so-called L s- 
channel QQ- cut’ was considered in mnanu. For a more detailed review see, for example, df 

I 11 the present paper we discuss all the different possibilities of colour-singlet y c (l + ) and 
yfc(l + ) inclusive production. We will not discuss here the colour-octet contribution, in particu¬ 
lar, since at the moment the importance of colour-octet contributions to the y c (l, 2) production 
remains open, see for instance m HU HU HE] . Note that Ref. [HU presents the complete NLO 
NRQCD predictions for the polarized y c (l,2) production at medium and high p t , which when 
compared to forthcoming LHC measurements could allow a detailed probe of the validity of 
NRQCD and colour-octet mechanism. 

We consider the lowest a s order process and, moreover, if in some kinematical domain we 
find that the cross section is enhanced by a large logarithm (either of virtuality, lng 2 , or of 
energy, Ins), then we will focus on this leading logarithm (LL) contribution. We emphasize 

Experimentally it is easier to measure such ratios since various theoretical and experimental uncertainties 
cancel out. As far as we aware, currently there are no high-energy data at low p t and no data on the rapidity 
distribution da/dy or on the total cross sections for y c ,b(l) except for the relatively low energy data in 7r Be 
collisions. Note that the p t of the J/%1) and not that of y c is measured, but due to the large J/if) mass (close to 
the Xc mass) the p t of J/ip is rather close to that of y c , and qualitatively reproduces well the behaviour with 
respect to the pt of y c 

3 Note that there is an indication for a rise of the ratio i? 2 i at low pt in the LHCb measurement of y c at 7 
TeV (B]. Therefore it is desirable that the LHCb collaboration perform a new measurement of this ratio at low 
p t during run-II of the LHC. 
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that the aim of the paper is not to present a precise quantitative prediction^} but rather to 
compare the role of different mechanisms of y(l + ) production in high energy hadron-hadron 
collisions. 

We note that Xb production has advantages over y c both from the theoretical and experi¬ 
mental viewpoints. Due to the larger mass of \b (i) the perturbative QCD approach is better 
justified and (ii) it is easier to detect the more energetic decay products. 

In Sect. 2 we first consider the cross section caused by the gg —> y(l + ) + g subprocess, and 
secondly due to the different virtualities of two incoming gluons. Then in Sect. 3 we discuss the 
production via three-gluon fusion where a pair of f-channel gluons represents the reggeization 
of one incoming gluon. Next, in Sect. 4 we study a most interesting possibility when the two 
gluons come from two different parton cascades. This mechanism for x(l + ) production has not 
been studied before. At asymptotically high energy this should be the dominating contribution. 
An interesting fact is that in such a case the cross section in the central region is expected to be 
smaller than that near the proton fragmentation domain. We give our conclusions in Sect. 5. 


2 Lowest order %(1 + ) production 

In this section we study production by the two-gluon initiated states shown by ‘1 and 2’ in Fig. [l] 
We show separately the contribution where all three gluons couple to the heavy quark loop 
(Fig.la), and the contribution where only two f-channel gluons couple to the heavy quark loop 
(Fig.lb) which may not vanish in for x(l + ) if the virtualities of these two gluons are different 
(Fig.lb). This contribution may be considered separately and formally called ‘production’ via 
two gluon fusion”. However actually this is just part of the whole gg —> y(l + ) +g cross section. 

Another point is that it is natural to consider for this part the role of f-channel gluon 
reggeization, which is the simplest example of x(l + ) production via 3-gluon fusion g + (gg) —>■ 
X(l + ). For this reason we will describe the contribution of Fig. lb in more detail in Section 
and the reggeization in Section [3} 

Finally, in Fig.lc we show the process initiated by a highly virtual s-channel gluon. 

2.1 The gg —>• x(l + ) + 9 process 

The simplest possibility to overcome the Landau-Yang selection rule is to create the x(l + ) 
meson together with an additional gluon. The corresponding ‘hard’ cross section was calculated 
in [TO] 

d°(gg -f x(l + ) + g) _ 12vr a 3 s R' 2 
dt ~ M 3 s 2 

P 2 [M 2 P 2 (M A - 4P) + 2 Q(-M 8 + 5M 4 P + P 2 ) - 15 M 2 Q 2 } 

(Q - M 2 P) 4 ’ ^ 

4 Precise results would require dedicated higher-order calculations. 
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Figure 1: Subprocesses for X(l + ) production: (a) the standard gg —> x(l + ) + 9 process where 
three on-mass-shell gluons couple to the heavy quark loop; (b,c) where one virtual gluon and one 
on-mass-shell gluon couple to the heavy quark loop — the virtual gluon has q 2 < 0 for (b) and 
q 2 > 0 for (c). 


where s, t, u are the Mandelstam variables for the hard subprocess. The quantities P = st + 
su + ut and Q = stu ; and M is the Xe(l + ) or X&(1 + ) meson mass. 

In eq.Q the nonrelativistic wave function of the meson is assumed. The derivative of the 
P- wave wave function at origin is denoted as R'. In our calculations for x c (l + ) we use the values 
R' 2 /M 2 = 0.006 GeV 3 (as in (103) an d the QCD coupling a s = 0.335 which provide reasonable 
widths of the x c (0,1, 2) mesons calculated accounting for the a s corrections (see e.g. [T7]). For 
Xfc(l + ) we take R' 2 = 1 GeV 5 similar to that in [7], which is consistent with the potential model 
results, and a s = 0.22 corresponding to the higher scale appropriate for \b production. 

In comparison with the natural parity X c (0 + ) and x c (2 + ) mesons, which can be formed via 
on-mass-shell gluon-gluon fusion with the cross section a on a 2 , the cross section ([Tj) contains 
and extra power of QCD coupling (not accompanied by any large logarithm). Therefore gg —)• 
x(l + ) + g can be considered as a NLO production process. 

Expression (Jl]) has to be convoluted with the incoming parton distributions and integrated 
over t and s. Since we are looking for the inclusive x(l + ) production we fix the rapidity, Y, 
of x(l + ) meson and choose to integrate over the variables y and pt, where pt is the transverse 
momentum of the final gluon while y is the rapidity separation between this gluon and the 
x(l + ) meson. It is easy to check that the corresponding Jacobian J — 1 (see e.g. [TEj)- Thus 
the cross section of x(l + ) production in proton-proton collisions reads 

= Jx\g{x\, Vf) ^" ^ + x 2 g(x 2 , y F )dydp 2 , ( 2 ) 

where g(x,/j,F ) is the density of gluons which carry a fraction x of the momentum of the 
incoming proton (measured at factorization scale fip). The variables 

aJi ,2 = {m t + p t e ±y )e ±Y / (3) 
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where m 2 = M 2 + p 2 . The hard cross section da/dt is given by Q with 

s = ml + p 2 t + r m t p t (e y + e ~ y ) (4) 

t = -p 2 t - m t p t e v (5) 

u = -pi - m t p t e~ y . (6) 

The expected cross section is shown in Table 1 where the LO MSTW2008 [19] PDFs were used. 


2.2 Production via g* + g —>■ x(l + ) fusion 


Here we consider the possibility of y(l + ) production via ‘virtual+reah gluon fusion processed 
which are shown in Fig. [TJ These processes are not forbidden thanks to the different mass/virtualities 
of the initial gluons. The corresponding ‘hard’ cross section can be extracted from ([!]) using 
the ‘equivalent photon/gluon’ approximation [2T]. Indeed, in the limit of s M 2 , |t|, the 
dominant contribution comes from a diagram where first the incoming gluon with momentum, 
say, pi emits the final gluon p 3 and then the virtual gluon, g*, with momentum q = pi — p 3 
interacts with another incoming (quasi-real) gluon p 2 to produce the y(l + ) meson |^] 

Thus in the large s limit we may write the cross section ([!]) as the product of the virtual 
gluon flux, dN = (asNc/^^dz/z)dq 2 / q 2 times the elementary g* + g —> x(l + ) cross section. 
That is 


da 

dt 


dN 

= +9 ->• Xc(l + )) 


(7) 


s>M 2 ,|t| dq 2 

(The factor dz/z in dN corresponds to the integration over the rapidity separation y and is 
omitted here.) In this way we get 


v(9* + 9~> X(l + )) 


4vr 2 alR' 2 AM 2 — 2t 
-\t\- 


M 3 


mt 


( 8 ) 


where the virtuality, q 2 , of off-mass-shell gluon g* plays the role of t = q 2 in ([!]). 

Since the elementary cross section (J8]) vanishes as q 2 —» 0 we cannot consider the incom¬ 
ing gluon, with momentum q, as an on-mass-shell parton. The inclusive cross section should 
therefore be written in terms of the unintegrated gluon density which is defined in such a way 


that 


xg(x,g F ) 



fg(x, f.VF) 



(9) 


Thus we obtain 


da 

dY 


J x x g(x i, g F )d-(g* + g x(l + ))f g {x 2 , q 2 , g F ) 



+ (x 2 Xi) , 


( 10 ) 


5 The production via g* + g fusion was discussed in detail in [20. 

6 Recall that in terms of the unintegrated gluon density fg(x,qt,y F ) the value of xg(x 1 y F ) is given by the 
logarithmic integral xg(x,y F ) = f yF f g (x,q t ,y, F ) dq 2 /q 2 . That is, each function xg contains a large logarithm. 
In this logarithmic integration the virtuality of the initial gluon is small, q 2 ~ q 2 <C yip. So this gluon may be 
considered as an on-mass-shell particle. 


5 






X 



&G - 


Figure 2: y(l + ) production via the fusion of an incoming gluon with a pair of gluons which form a 
Regge trajectory a G - 


with ,x' 1)2 = e ±y y/ (M 2 — q 2 )/s. 

Note that the dq 2 /q 2 integral does not now have a logarithmic structure since the ‘hard’ 
cross section ([8]) contains a factor t = q 2 . In other words expression (10) should be considered 
as a NLO contribution (in comparison with LO x(0 + , 2 + ) meson production) at the same level 
as the cross section (J2]) . Moreover, it should not be considered as a new contribution - it is just 
a part of the whole cross section (J2|. 

The value of contribution (10) shown in Table 1 is obtained by calculating the unintegrated 
gluon density starting from the MSTW2008 [19] PDFs based on the KMR/MRW last-step 
prescription [22j [23]. In this prescription we have used the LO splitting functions but kept the 
more exact NLO kinematics. As it was shown in [23j this provides an accuracy close to that 
given by the NLO prescription. Note that here we have used the gluons unintegrated over the 
virtuality, q 2 and not over the transverse momentum squared, q 2 (see [23J for details). 

The contributions Q and (10) were evaulated using the MSTW2008 LO PDFs for the 
integrated gluons xg. For the very low q 2 < Ql we assume the saturation-like behaviour 
xg{x, q 2 ) = xg(x , Q 2 0 )(q 2 /Q 2 0 ). 


3 Production via gluon reggeization 

An attractive possibility to organize the fusion of one gluon with the pair of incoming gluons, 
that is the g + (gg) —> x(l + ) subprocess, is to consider the contribution coming from gluon 
reggeization, see Fig. [2} On one hand the gluon trajectory is described by the diagrams where 
the initial t-channel gluon is replaced by the exchange of two f-channel gluons. On the other 
hand this contribution is enhanced by large logarithms of the proton-proton energy, or to be 
more precise - by a 1/ooq factor, where loq denotes the shift of the position of the BFKL vacuum 


6 






singularity «bfkl = 1 + u; 0 . Recall that thanks to the bootstrap condition [23j if we account for 
the gluon trajectory then we include all the contributions (of the antisymmmetric colour-octet 
states, which are of interest here) enhanced by the 1/ujq factor, that is - by the large leading 
logarithm of the energy. 

Let us explain in more detail how the gluon trajectory is built up. Consider the diagrams 
where the additional ^-channel gluon is added to the usual ladder diagram (inside the same 
parton cascade). The simplest such diagram (in Feynman gauge) is where the additional gluon 
is between the y c and the nearest s-channel gluon of the ladder (that is, between the c-quark 
loop and the gluon 3 in Fig.lb). The contribution of this diagram is enhanced by two logarithms. 
One logarithm, lng 2 , comes from the integration over the k t of the new f-charmel gluon. The 
other logarithm comes from the integration over the longitudinal component of the new gluon’s 
momentum, which corresponds to the integration over the mass of the intermediate s-channel 
gluon which emits the new t-channel gluon. This longitudinal logarithm is actually equal to 
the rapidity separation between the y c meson and the nearest s-channel gluon. This rapidity 
separation is driven by the intercept (the ^-dependence) of the parton cascade. In the BFKL 
approach this logarithm is equal to 1/ujq. Moreover, in terms of the BFKL amplitude, the 
coherent sum of such diagrams (where the new f-channel gluon couples to the different s- 
channel gluons in the ladder) is described by the gluon trajectory 

«c(g 2 ) = • (11) 


Note that in this approach (just as in two gluon fusion g* + g —> x c (l + )) the corresponding 
contribution vanishes when the ‘reggeized’ gluon virtuality q 2 —» 0. Indeed, due to gauge 
invariance, the whole set of diagrams which describe gluon reggeization can be reduced to a 
one-gluon loop inserted in the place of the original gluon propagator (see, for example, [23]). 


Thus the result looks like cross section (10) multiplied by the gluon trajectory ac(q 2 ) and by 
the 1/uo factor. 

Moreover, contrary to inclusive J/if production ( where the gluon pair should be in a 
symmetric colour-octet state and the analogous contribution is imaginary), in the case of \ / (l + ) 
we deal with the antisymmetric colour octet — that is, with a true gluon trajectory of negative 
signature. Therefore the logarithmically enhanced contribution is real and interferes with the 
lowest a s order g* + g — > Xe(l + ) amplitude considered in the previous section. 

Thus, finally, to lowest order in a s , the cross section (10) should be multiplied by a ‘double 
logarithmic’ factor A, where 

. a s N r 


A = 


1 + 2 


2tUjJc 


■ ln («l/«2 


( 12 ) 


Here the q 2 logarithm comes from the integration over the transverse momentum k t in the loop 
corresponding to the gluon trajectory. In the region of k 2 -C q 2 , the integral takes the form 


dk 2 


(13) 
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Note that the amplitude for unnatural parity 1 + meson production 


9(Qi) +9(42) -► X(l + ) 


( 14 ) 


is antisymmetric with respect to the permutation of the two gluons. Therefore in the amplitude 
(14) we have the logarithm of the ratio q\/q This result solves the problem of the infrared 
cutoff at low k t in the integral (13). In the pure symmetric configuration (with q\ = q%) we get 
zero. This fact in some sense is similar to the Landau-Yang selection rule - a spin 1 + particle 
cannot decay into two identical transverse gluons [[} 


Recall that in the case of y(l + ) production the second f-channel gluon (whose distribution is 
written here in terms of the ‘gluon reggeization’ trajectory (11)) is not a virtual loop correction 
to the main amplitude. Rather, this is a particular new channel to produce a spin 1 + meson. 
For this reason it should not be combined with the real, s-channel gluon emission. The infrared 
divergence in (12) is cancelled between the diagrams with the upper (qq) and the lower (q 2 ) 
gluon insertions; and not between the virtual loop (reggeization) and the real gluon emission^ 


The results of the numerical estimate of the corresponding order of a s contribution^] are 
presented in the fourth column of Table 1. Here we have used the value of loq = 1/4 which 
is close to that expected for the BFKL pomeron after the resummation of the next-to-leading 
log corrections [26]. Recall that at the lowest a s order, the ‘reggeized induced’ contribution is 
proportional to a s /ojQ. We choose reasonable values of cuo and a 8 to indicate the possible size 
of the effect. 


As expected, the result is negative since at q 2 < 0 the gluon trajectory is shifted to lower 
values of a G (q 2 ) < 1. However due to the cancellation between the qq and q 2 terms the whole 
contribution is not too large inspite of the 1/loq enhancement. 

Recall however that at the present stage we account for the order of a s reggeized correction 
only. When the ’’correction” becomes large the higher order a s terms become important re¬ 
placing effectively the first and negative a s contribution (let denote it as — Sr) by the positive 
exponential factor like exp(— Sr). 


4 Production via two parton cascades 

Here we consider the situation when all three (quasi-real) gluons couple to the heavy quark loop 
directly. The most interesting possibility is to form the incoming (gg) pair taking the two gluons 

7 Recall that as we are looking for the leading logarithm the gluons may be considered as quasi-real, transverse 
particles. 

8 Just as in the case of the colourles Higgs boson production, the infrared divergence caused by real emission 
is cancelled by the true reggeization diagram where (in Feynman gauge) the additional f-channel gluon couples 
to the s-channel gluons above and below the colourless boson. 

9 The contribution of the terms oc In qf and oc In are calculated separately. In the first case the gluon qi 
is written in terms of the unintegrated distribution, f g , while for gluon q 2 we may use the integrated gluon 
distribution x 2 g(x 2 , /j,f), and vice versa. 
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Figure 3: The diagram for the cross section for x(l + ) production via the fusion of two gluons from 
different parton cascades, where pa and p B are the four momenta of the incoming protons (not 
shown). The diagram is for the cross section, AA* , so the particles intersected by the dashed line 
are on-mass-shell. 


from two different parton cascades, see Fig. [3j The probability to find the corresponding pair 
is given by the product of the gluon densities, x 1 g(x 1 ) x X 2 g(x 2 ), multiplied by the probability 
that the two cascades overlap in transverse space. Since in the low x region (which is relevant 
at high energies) the parton density grows as a power of (1/x) - that is xg(x) oc x~ x - this 
contribution will dominate asymptotically at very high energies. 

In Double Parton Scattering (DPS) the probability of cascade overlap in transverse space 
is given by a factor 1/<7^ S (see e.g. [27j), where 


DPS 


cr. 


eff 


- / . 


(15) 


Actually this factor results from the integration (15) over the (‘pomeron’) loop formed by the 
two gluon cascades. The integral is driven by the proton form factor F(t), that is by the t ~ q 2 
dependence of the proton-pomeron vertex. However, while for DPS the loop contains four form 
factors (leading to F 4 in (15)) in our case we have only two form factors (we need two cascades 


on one 


side of Fig. [3] only). Therefore (assuming exponential dependence) we have 


cr. 


x _ 


eff 


1 

= -cr. 


DPS 

eff 


(16) 


We emphasize that the typical value of \t\ in the pomeron loop integration simultaneously plays 
the role of the lower limit for the factorization scale, p? F . In terms of <J^ PS (and assuming the 
exponential f-behaviour) this limit is p 2 F > (|t|) = 47 t/cj^ ps . 

Let us denote the momenta of the three gluons as F (i — 1, 2, 3). In terms of the incoming 
proton’s momenta, pa and p B we may write 


fci ~ ap A + k lt 
h ~ PiPB + ku for i = 2,3. 
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(17) 

(18) 








To calculate the corresponding matrix element we use the gauge invariance condition (Ai^k^ = 
0) and replace the proton’s momenta Pa,Pb transferred through the spin part (or numerator) 
of gluon’s propagators by fat/a and k lt / fa, respectively. Note that in order to keep the leading 
(DGLAP) logarithms in k t integrals we retain the lowest power of k lt in the matrix element p°j 
Thus the matrix element for ggg —> y(l + ), corresponding to the heavy quark loop in Fig. [3j 
takes a rather simple form 

M = -iBf abc ^~, (19) 

N C M v ’ 

where f abc is the antisymmetric colour tensor, N c — 3 is the number of colours, and g is the 

QCD coupling (a s = g 2 / 47t). The basic amplitude A s contains a spin part (given by the trace 

around the quark loop) and a part corresponding to the propagator poles. It is of the form 


where the trace gives 


and the poles are 


A s — 16(aiZi + CL 2 Z 2 ) 

(20) 

ai = (e a • k-st) e Q/ 3 7< 5 k% t p p A pi ej 

(21) 

02 = (e a • k 2t ) e a/3l 8 k% pP A pi e s x , 

(22) 

_ 2 a(3 2 3 s 

a(2z - 1)^235 - \k u - fat | 2 - M 2 

(23) 

2a/323 s 

a/3 23 s(l - 2z) -\fa t - fat\ 2 ~ M 2 ’ 

(24) 


where the denominators include the contribution of the longitudinal (afa^sfa — 2 z)) and trans¬ 
verse (—\fat ~ fat] 2 ) components of the square of the momentum. Here e x and e a are the x(l + ) 
and A;i-gluon polarization vectors. Also we have introduced the relative momentum fraction 
z = fa/(fa + fa) where 


@23 — @2 + @3 — (AT 2 + \kit + fat\ 2 )/oiS. 


(25) 


The value of /? 2 3 is fixed by the y(l + ) meson mass and its rapidity. Finally, the normalization 
constant B in (19) is given by 


B = 


3 R' 2 
7 rM 3 


(26) 


Recall that in the diagram for the cross section, Fig. [3j we have four /-channel gluons from 
the proton ps side [^] These 4 gluons form three loops of integration. The integrals over 


10 Together with the denominators (1 /kf) of the gluon propagators this will give in the cross section the 
logarithm j k 2 t d 2 k lt /kf. 

11 In general, there may be a diagram where the gluon pair ( gg ) goes in the pb direction in one amplitude, A, 
but goes in the other, pa, direction in the amplitude A*. Such a contribution either corresponds to the three 
gluon singularity (in the u plane) which has a lower intercept (wo ~ 0 or even negative - this contribution is 
small): or, dealing with the BFKL pomeron, we have to consider this (gg) pair as gluon reggeization which (in 
the lowest a s order) was already considered in Sect.3. 
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the transverse momenta components have already discussed: two loops gives the logarithmic 
integrals f k 2 t d 2 k it /kf corresponding to the gluon PDF given by the parton cascade i, while 
the third integral (over t) is limited by the proton form factor. The longitudinal momentum 
component in the central loop is fixed by the condition that the y(l + ) meson (and other 
secondary partons) should be on-mass-shell. Finally, we have the loop integrations over the 
gluon longitudinal momentum in each amplitude, A and A*. These integrals can be written in 
terms of z and can be closed on the quark pole[^] that is, in our non-relativistic approximation 
(for the y(l + ) wave function) we get an intermediate state with both quarks on-mass-shell. 
The pole position is 


2{k 2 t • &3t) 

M 2 


(27) 


So, finally, the matrix element (19) reads 


M 


32nBf abc 


g 3 (a 1 + Qg) 

N C M 


(28) 


To obtain the cross section we have to sum over the y(l + ) polarizations, e x , average over 
e a and integrate over the gluon transverse momenta. A problem is that according to (27) the 
value of z depends on these momenta and we cannot use the usual, integrated PDFs; each value 
of kn corresponds to its own z. Therefore we have to use the unintegrated distribution in the 
same way as described in Sect. 2.2. Moreover in our calculation we have to keep only the poles 
with the positive z < 1. Otherwise one of the gluons (z or 1 — z) will have negative energy 
and so must be considered as an outgoing and not an incoming particle. Such a configuration 
describes the gg —> y(l + ) + g subprocess which has already been discussed in Sect. 2. 


Thus the ‘DPS’ component of inclusive y(l + ) cross section (which we call u DPS ) reads 


da DPS 

dY 


n 3 a 3 s R' 2 32 
uD ps M 9 ~9~ 



dk lt 

k 2 k 2 
^2 1 'D t 


otg(a, M 2 ) 


J tjc ^ 1 + C ° s2 ^^ X2 = ^ 23 ( 1 ~ k %t)fg( x 3 = ZP 23 , kl t )Q(z)Q(l - z) + (a ++ /3 23 ) , (29) 

where z is given by ( [27] ) and where we take the factorization scale Hf — M. Since, within the 
Leading Logarithm approximation, both k\ v k^ t M 2 , we may replace in the first argument 
of f g ((3 23(1 — z ),...) the value of (1 — z) by 1, that is / ff (/? 23 (l — z), ...) becomes f g {/3 2 3, ...). The 
origin of cos 2 (j) in the cross section comes from the convolution of a\ and a 2 . 

The predicted ‘DPS’ part of the x(l + ) production cross section is presented in the fifth 
column of Table 1. In the numerical calculation we have used MSTW2008LO partons [T9] and 
(x^ PS = 10 mb [29] . 

12 In general, in our case with a L = 1 wave function, some part of the contribution may contain poles of 
second-order (see e.g. [IS]). However, as far as we keep just the leading logarithms in k t , that is the lowest 
power of k t , the residue of these second-order poles vanishes at z = 0 or 2 = 1. Therefore actually we deal only 
with simple first-order poles. 
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Y 

gg ^x + g 

:k . 

g g^x 

G-regge 

DPS 

0.0 

9.8 (17.2) 

6.3 (10.5) 

0.0 (0.0) 

1.1 (2.2) 

1.0 

9.5 (16.8) 

6.2 (10.3) 

-1.0 (-1.3) 

1.2 (2.4) 

2.0 

8.7 (15.6) 

5.6 (9.6) 

-1.9 (-2.5) 

1.4 (2.9) 

3.0 

7.4 (13.7) 

4.7 (8.5) 

-2.5 (-3.3) 

1.7 (3.6) 

4.0 

5.6 (11.3) 

3.6 (7.1) 

-2.9 (-3.8) 

2.1 (4.7) 

5.0 

3.6 (8.5) 

2.1 (5.4) 

-2.9 (-4.1) 

2.3 (5.9) 

6.0 

1.5 (5.1) 

0.6 (2.9) 

-2.4 (-4.5) 

1.8 (4.6) 


Table 1: The cross section da/dY in fib for producing Xc(l + ) by the various mechanisms at 7 (13) 
TeV. Note that the g*g —* % is already included in the gg —>■ y + g contribution, and has only been 
considered separately to facilitate the study of the Regge contribution. 

5 Discussion 

As seen from Tables [Tj and [2j in the central rapidity region (Y = 0) the main contribution to 
y(l + ) production comes from the most trivial gg —> / \(1 + ) + g subprocess, even at 13 TeV. 
This contribution decreases with |Y|, that is, towards the edges of the available rapidity space. 

The expected correction due to gluon reggeization vanishes at Y = 0, in accord with the 
generalized Landau-Yang selection rule, which is valid for the symmetric configuration. More¬ 
over, for Y yf 0 it is negative, and at large rapidity Y it becomes quite large in comparison with 
the original (without reggeization) contribution. The absolute value of the Regge contribution 
increases with |Y| (up to Y ~ 5) making the overall rapidity (Y) distribution a bit narrower. 
We emphasize, however that the results of Tables [l] and [2] are presented just for illustration 
purposes, and include the 0(a s ) reggeized correction only. When the ‘correction’ increases, 
the higher-order a s terms become important replacing effectively the first and negative a s 
contribution (which we denote as — 5r) by a positive exponential factor like exp(— Sr). 

The DPS contribution, which originates from the fusion of gluons from two parton cascades, 
reveals quite a different behaviour. Due to the growth of the low x gluon density this part of 
cross section increases with the initial energy and/or with |Y|. The DPS mechanism for y c (l + ) 
production is seen to provide about 30% of the cross section for Y = 4 to 5. At very high 
energies (asymptotically) this contribution starts to dominate. 

The values of the Xb(l + ) cross sections are much smaller, but due to the larger \b mass the 
perturbative predictions are better justified. Note also the strong correction caused by gluon 
reggeization, which arises since the virtuality of the gluon is larger. 
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Y 

99 X + 9 

jk . 

9 9 

G-regge 

DPS 

0.0 

37 

30 

0 

2.2 

1.0 

36 

29 

-10 

2.4 

2.0 

31 

25 

-18 

2.7 

3.0 

24 

18 

-22 

3.1 

4.0 

16 

11 

-23 

3.1 

5.0 

8 

3 

-18 

2.4 


Table 2: The cross section da/dY in nb for producing Xb(l + ) by the various mechanisms at 13 
TeV. Note that the g*g —> \ ' s already included in the gg —> x + 9 contribution, and has only been 
considered separately to facilitate the study of the Regge contribution. 
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